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We use Minkowski Functionals to explore the presence of non-Gaussian signatures in simulated 
cosmic microwave background (CMB) maps. Precisely, we analyse the non-Gaussianities produced 
from the angular power spectra emerging from a class of inflationary models with a primordial step¬ 
like potential. This class of models are able to perform the best-fit of the low-^ ‘features’, revealed 
first in the GMB angular power spectrum by the WMAP experiment and then confirmed by the 
Planck collaboration maps. Indeed, such models generate oscillatory features in the primordial 
power spectrum of scalar perturbations, that are then imprinted in the large scales of the GMB 
field. Interestingly, we discover Gaussian deviations in the GMB maps simulated from the power 
spectra produced by these models, as compared with Gaussian AGDM maps. Moreover, we also 
show that the kind and level of the non-Gaussianities produced in these simulated GMB maps are 
compatible with that found in the four foreground-cleaned Planck maps. Our results indicate that 
inflationary models with a step-like potential are not only able to improve the best-fit respect to the 
AGDM model accounting well for the ‘features’ observed in the GMB angular power spectrum, but 
also suggesting a possible origin for certain non-Gaussian signatures observed in the Planck data. 

PACS numbers: pacs 


I. INTRODUCTION 

The Standard Cosmological model is efficiently sup¬ 
ported by the inflationary paradigm to explain the flat¬ 
ness and homogeneity of the observed Universe. At the 
same time, the inflation model provides also an elegant 
mechanism to produce the primordial curvature pertur¬ 
bations that generate the seeds for the formation of struc¬ 
tures. 

The most recent cosmic microwave background (CMB) 
data by the Planck satellite [T] are in excellent agree¬ 
ment with the assumption of adiabatic primordial scalar 
perturbation with nearly scale-invariant power spectrum, 
described by a simple power law with spectral index rig 
very close to (albeit different from) unity [2] . It would be 
produced in the simplest inflationary scenario, in which 
a single, minimally-coupled scalar field slowly rolls down 
a smooth potential. 

In spite of this, models that accounts for the localised 
‘features’ in the primordial power spectrum provide a 
better fit to the data with respect to a smooth power-law 
spectrum [Ml]. Features in the primordial power spec¬ 
trum can be generated following departures from slow 
roll, that can happen in more general inflationary mod¬ 
els with a symmetry-breaking phase transition. It is the 
case of the inflationary models with a parametrized step 
in the primordial potential, which show oscillations in 
the power spectrum of curvature perturbations, localised 
around the scale that is crossing the horizon at the time 
the phase transition occurred [IMS]. 

The main tools for analyse (and constrain) the infla¬ 
tionary models are the current CMB temperature power 
spectrum and bispectrum [TBUT^ . Another important 
approach to study viable inflation models is through 
the non-Gaussian signatures produced during the infla¬ 


tionary phase, that left imprints in the CMB tempera¬ 
ture fluctuations. Recent constraints on primordial non- 
Gaussianities (NGs) found in the Planck Collaboration 
analyses |20l [21] severely constrain those of local type 
(hereafter local-NG), supporting the simplest inflation¬ 
ary model based on a single minimally-coupled scalar 
field. However, such analyses do not rule out inflationary 
models producing other types of primordial NG, when¬ 
ever consistent with such constraints. In this scenario, 
inflationary models with a step-like feature in the infla- 
ton potential deserve particular consideration [22] [23] . 

Non-Gaussian signals, primordial or not, appear mixed 
in the CMB Planck maps, with each phenomenon con¬ 
tributing with their own signature. Various statistical 
methods are proposed in the literature [20] M [HHH] 
in order to detect all the potential NG components in 
the CMB, to measure their intensity and their angular 
scale dependence. Here we use the Minkowski Function¬ 
als (MFs) [30] as statistical estimators to look for non- 
Gaussian features in several sets of simulated maps. Be¬ 
tween others, we analyse Monte Carlo CMB maps, seeded 
by the CMB angular power spectrum originated from in¬ 
flationary models with a step-like feature in the inflaton 
potential, and compare them with simulated CMB maps 
based on the ACDM concordance model, for which we 
use the current data of the Planck collaboration m- 

This paper is organised as follows: in Section [TT| we 
briefly recall the inflationary formalism and introduce the 
class of models with a step-like potential. In Section [FH] 
we briefly introduce the MFs as statistical estimators of 
NGs. In Section [IV] we present the results of our analy¬ 
ses, and finally, in Section |V| we display our concluding 
remarks. 
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II. INFLATIONARY PERTURBATIONS 

In this section we first introduce the formalism to cal¬ 
culate the power spectrum and bispectrum from a given 
inflationary potential V{(j)). Then, we describe the de¬ 
partures from the slow-roll approximation and introduce 
the inflationary class of models with step. 

The Lagrangian for the inflaton scalar field 0, mini¬ 
mally coupled to gravity is 

C^ = \g^''d^4>d,4>-V{4>), 

with V{(j)) its potential. The equation of motion for a 
homogeneous mode of the field (l){t,x) is described by 
the Klein-Gordon equation 


0 + 0. (1) 

This is a familiar equation for a free scalar field with an 
extra term (3i70, with H the Hubble factor) that comes 
from the use of the Friedmann-Robert son-Walker (FRW) 
metric in the Lagrangian. The Klein-Gordon equation 
and the first Friedmann equation 
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where mpi is the Planck mass, are used to determine the 
background dynamics for both the Hubble parameter and 
the (unperturbed) inflaton field (j). 

In order to study the evolution of the primordial scalar 
fluctuations, we introduce the comoving curvature pertur¬ 
bation IZ = and the two-point correlation function 


(7^fc7^fe.) = i2TTf5ik + k')VTt{k). (3) 

with Vn the curvature power spectrum. If IZ is Gaus¬ 
sian, the Eq. depends only on the magnitude of k if 
we use the isotropy assumption. If we also assume the 
homogeneity, we get the relation {TZk) = 0. 

Using the gauge-invariant quantity u = a5(j)^ where a 
is the scale factor, we can define 


u = —zlZ , 

with 2 ; = a^/H. The Fourier components of u obey to 
the equation 

u'k + Wfe = 0, (4) 

with a prime denoting a derivative with respect to con¬ 
formal time T]. 

The z" jz evolution depends on both the dynamics of 
the Hubble parameter and the unperturbed inflaton field, 
governed by the Friedmann equation (Eq. (§) and the 
Klein-Gordon equation for 0 (Eq. §). Integrating the 
Eq. 0- we get the U]z[r]) for free-field initial conditions. 


The power spectrum of the curvature perturbations 
Vp is related to u and 2 ; through 


VTt{k) 


27r^ 2 ; 


( 5 ) 


evaluated when the mode crosses the horizon. Assuming 
gaussianity and adiabaticity, this quantity contains all 
the necessary information for a complete statistical de¬ 
scription of the fluctuations. Instead, for non-Gaussian 
fluctuations the higher-order correlations contain addi¬ 
tional information. These require more statistics and are 
therefore more difficult to measure, especially at large an¬ 
gular scales where cosmic variance errors are significant. 
In case of NG, the fluctuations have a non-zero three- 
point correlation function Bn{ki; k 2 ; ks)^ named bispec¬ 
trum. One of the simplest form of local-NG is described 
by the parameterization [33]: 


nix) = Ugix) + 5/nl [ngix)^ - ingix)^)] , (6) 

where IZg is Gaussian. In this model for /ocaZ-NG, the 
information is reduced to a single number /nl- 

The level of NG generated during the inflation de¬ 
pends from the inflationary model considered. In case 
one assumes a inflationary model with a single scalar 
field, canonical kinetic terms, initial adiabatic (Bunch- 
Davies) assumption, and slow-roll condition, then a tiny 
level of NG is predicted (see, e.g., [23|). Thus, a robust 
detection of |/nl| > 1 would rule out this class of models, 
and vice versa. 

Indeed, when one of these assumptions is violated, a 
large amount of NG is expected |32|. In this case, such 
models can generate a detectable and unique signal of NG 
of local type, where the bispectrum amplitude is maxi¬ 
mized for the condition /ci ^ /c 2 ^ called squeezed 
triangle configuration. 

In this scenario, inflationary models with a step-like 
feature in the inflaton potential would have produced de¬ 
partures from Gaussianity during the inflationary phase, 
leaving detectable signatures in the CMB data ESEa], 
as we will examine. 


A. Inflationary step-like models 

In order to solve the Eqs. 0-0 is useful to intro¬ 
duce the slow roll approximation, which assumes a slowly 
varying inflaton field with (1/2)0^ ^ ^(0)- fhis way 
we can neglect the kinetic term in the Friedmann equa¬ 
tion and the acceleration term in the Klein-Gordon equa¬ 
tion. 

In some inflationary models the slow roll assumption 
could be relaxed for a few time, assuming the system 
starts in a state where the slow roll conditions are fulfilled 
(in order to give it enough time to reach the inflationary 
attractor solution), and the system returns to the slow 
roll regime at a later time. 
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The interruption of slow-roll leave possible detectable 
traces in the primordial power spectrum. Specifically, 
wavelengths crossing the horizon during this fast-roll 
phase will be affected, leading to a deviation from the 
usual power-law behaviour at these scales. 

The brief violation of the slow-roll condition can be 
modelled, in a single-field model, by adding a local fea¬ 
ture, such as a step m to an otherwise flat infiaton 
potential. This step can be regarded an effective field 
theory description of a phase transition in more realistic 
multi-field models [34], which may arise naturally in, e.g., 
supergravity- [35| or M-theory-inspired inflation models 

ES]- 

In this work we choose to analyse models with step-like 
features in the inflationary potential using the formalism 
by Adams et. al na that add a step feature to a V (0) = 
j2 chaotic potential, i.e., by considering a potential 
of the form 




1 -h ctanh 



( 7 ) 


where b is the value of the field where the step is located, 
c is the height of the step and d its slope. 

The spectrum of primordial perturbations resulting 
from the potential in Eq. 0 can be calculated as out¬ 
lined in the previous section, and is found to be es¬ 
sentially a power-law with superimposed oscillations, as 
showed in Fig. 0. The oscillations are localized only 
in a limited range of wave-numbers (centred on a value 
that depends on h) so that asymptotically the spectrum 
recovers the familiar form typical of slow-roll in¬ 

flationary models. 

The signatures of this class of models on the CMB tem¬ 
perature power spectrum and bispectrum HHni na IS3- 
m] and the tensor spectrum [42||43] were studied in detail 
using the current data, showing that these models im¬ 
prove the values with respect to the ACDM model. 
A similarly featured potential was also implemented in 
other context (see, e.g., j44l[46] b Here we go beyond 
these results, looking at the non-Gaussian signatures of 
these inflationary models imprinted on the CMB maps. 


III. ANALYSIS METHOD 

We performed a set of analyses using four Minkowski 
Functionals in order to search for non-Gaussian contribu¬ 
tions potentially present in simulated and Planck CMB 
maps. We are primarily interested in answering the fol¬ 
lowing related questions: (i) are there non-Gaussian con¬ 
tributions left in simulated CMB maps generated by the 
angular power spectra of inflationary step-like models? 
If the answer is yes, our second question is: (ii) are 
these non-Gaussian signatures compatible with the last 
released foreground-cleaned Planck CMB maps? 

In fact, if one finds non-Gaussian contributions in sim¬ 
ulated maps produced from the CMB power spectra ob¬ 
tained from inflationary step-like models, then we can in¬ 



k(Mpc''') 


FIG. 1. Primordial power spectra for inflationary potentials 
with step. The red line shows the reference power-law spec¬ 
trum, while the blue line draws the best-fit spectrum for the 
inflationary step-like model with m = 7.5 x 10~®, using the 
best-fit values of 


fer that these NGs are of primordial origin. For this the 
pertinency of the second question, to discover whether 
these non-Gaussian features are present, or not, in the 
foreground-cleaned Planck maps. 


A. The Models 


We start our analysis from the results of the previ¬ 
ous work of Benetti m, using the following values for 
the cosmological and step parameters: 100 = 2.22, 

= 0.1212, r = 0.089, 100(9 = 1.0411, n, = 0.959, 
10^As = 2.20, b = 14.66, c = 1.4 x 10-^ d = 3.6 x 10-^ 
and we called this model as Model A. 

We also consider a second model, called Model using 
the same cosmological and step parameters values of the 
previous model with the exception of the step amplitude 
parameter c, that is four times higher than in the Model 

A. 


We performed a Monte Carlo Markov Chain analy¬ 
sis via the publicly available package CosmoMC m and 
used a modified version of the GAME ([48]) code in which 
we numerically solve Eqs. 0-0 using a Bulirsch-Stoer 
algorithm in order to theoretically calculate the initial 
perturbation spectrum Eq. 0 needed to compute the 
CMB anisotropies spectrum. The resulting CMB power 
spectra for the Model A (dashed line) and Model B (dot- 
dashed line) are showed in Figs. 0 and 0 , together with 
the featureless ACDM angular power spectrum (contin¬ 
uous line). 
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B. The Minkowski Functionals 

The Minkowski Functionals (MF) present many ad¬ 
vantages as compared with other statistical estimators, 
between them one can mention their versatility to de¬ 
tect diverse types of NG without a previous knowledge 
of their features, like their intensity or angular depen- 
dance, and also because they can be efficiently applied in 
sky patches or masked maps. 

All the morphological properties of a d-dimensional 
space can be described using d + 1 MFs m- In the case 
of a 2-dimensional CMB temperature field defined on the 
sphere, AT = AT(d, 0), with zero mean and variance <j^, 
this tool provides a test of non-Gaussian features by as¬ 
sessing the properties of connected regions in the map. 
Given a sky path V of the pixelized GMB sphere 5^, an 
excursion set of amplitude Ut is defined as the set of pixels 
in V where the temperature field exceeds the threshold 
that is, the set of pixels with coordinates (d, 0) G V 
such that AT(d,0)/cr = u > Uf. Each excursion set, or 
connected region, E, with u > and its boundary, 5E, 
can be defined as 

i: = {{e,4>)er\AT{e,4>)>iy<jo}, ( 8 ) 

as = {(0,0) e P I AT{e, 0) = jyao}. (9) 

In a two-dimensional case, for a region E C 5^ with am¬ 
plitude ut the partial MFs calculated just in Ri are: a^, 
the Area of the region described by E, the Perimeter, 
or contour length, 9E of this region, and the number 
of holes inside E. The global MFs are obtained calcu¬ 
lating these quantities for all the connected regions with 
height u > Uf. Then, the total Area A(z/), Perimeter L{u) 
and Genus G{n) are [2411491 l5Q] 

= J- ^ / Kdl = y^gi = Nhot - Neoid , (12) 

where d£L and dl are the elements of solid angle and line, 
respectively. In the Genus definition, the quantity n is 
the geodesic curvature (for details see, e.g., [50]). This 
last MF can also be calculated as the difference between 
the number of regions with v > Vt (number of hot spots, 
Nhot) and regions with n < Uf (number of cold spots, 

^cold)' 

The calculations of the MFs used here were done us¬ 
ing the algorithm developed by m and m- This code 
calculates four quantities, namely Vq = Vi = L{n)^ 

V 2 = G(z^), and V 3 = Nduster s{t^) which are the three 
usual MFs defined above plus an additional quantity 
called number of clusters, Nciusters{^)- The latter quan¬ 
tity is the number of connected regions with height u 
greater (or lower) than the threshold Uf if it is positive 
(or negative), i.e., the number of hot (or cold) spots of 
the map. 


C. Data description 

The five sets of simulated maps we investigate here 
were produced with Ngjq^ = 512, and in the case of the 
Planck GMB maps, they were degraded to this resolu¬ 
tion, using the HEALPix (Hierarchical Equal Area iso- 
Latitude Pixelization) pixelization grid [52]. In this sec¬ 
tion we give more details about these sets of synthetic 
maps as well as the Planck maps in analysis. 


1. The Monte Carlo CMB maps 

Eor our analyses we produce five sets of Monte Garlo 
(MG) GMB maps, random realisations considering a 
given GMB angular power spectrum. Each set is com¬ 
posed by 1000 GMB maps. The first three MG sets are 
seeded by: (i) the AGDM angular power spectrum, ac¬ 
cording to last Planck results (solid line in Eig. HI; (ii) 
Model A, corresponding to the best-fit of Planck angular 
power spectrum data using an inflationary model with 
step-like potential (dashed line in Eig. [^; (hi) Model 
B, which also uses an inflationary model with step-like 
potential but fits less accurately the Planck data (dot- 
dashed line in Eig. [^. 



Multfpole 1 

FIG. 2. Angular power spectra seeding the spectra of the sets 
of simulated CMB maps analysed here: the ACDM (solid blue 
line), the Model A (dashed red line; see text for details), and 
the Model B (dot-dashed orange line; see text for details). 

The last two sets of MG GMB maps are also generated 
using the AGDM power spectrum, but including different 
contributions of local-NG, namely, (iv) with /nl = 38, 
and (v) with /nl = 100, the last set. The first one is 
chosen to be consistent with recent constraints found by 
Planck [20] [21], that is, /nl = 38 ±18 (at 68% confidence 
level) for the large angular scales. Moreover, our choice 
foi* /nl = 100 aims to test also a local non-Gaussian con¬ 
tribution with higher amplitude in order to compare the 
signatures they produce in the MFs. These latter two 
sets of non-Gaussian GMB maps were generated from 
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a set of coefficients a£m derived from a combination of 
the multipole expansion coefficients {a^^} and 
corresponding to CMB Gaussian and non-Gaussian (of 
local type) maps, respectively. A set of each kind of 
these spherical harmonics coefficients, 1000 linear and 
1000 non-linear, were produced by [53], and are publicly 
availabl43 We combine these data as follows 

= apm +/nl , (13) 

which were then normalised by the Planck best fit power 
spectra m, rescaling by the ratio of the square root 
of the power spectra and directly by the ratio of the 
power spectra. This linear combination permits to derive 
GMB synthetic maps with an arbitrary level of of local- 
NG defined by any real value of /nl- 

As can be observed in Fig. the main differences be¬ 
tween the ACDM power spectrum and those from Models 
A and B are mainly concentrated at large scales, specif¬ 
ically for = 10 — 50. For this reason all our analysis 
will be done focusing on this range of the angular power 
spectra. That is, in all the MG GMB maps here anal¬ 
ysed we consider only the angular scales of this range of 
multipoles, i.e., ^ = 10 — 50 (see Fig. [^. 

2. The foreground-reduced Planck CMB data 

In February 2015, the Planck Gollaboration released 
the second set of products derived from the full Planck 
dataset Between them are the four GMB foreground- 
cleaned maps, namely Spectral Matching Independent 
Gomponent Analysis (SMICA), Needlet Internal Linear 
Gombination (NILC), Spectral Estimation Via Expecta¬ 
tion Maximization (SEVEM), and the Bayesian paramet¬ 
ric method Commander, names originated from the sep¬ 
aration method used to produce them [54]. The per¬ 
formance of the foreground-cleaning algorithms has been 
carefully investigated by the Planck collaboration with 
various sets of simulated maps. These high resolution 
maps have Ngjjg = 2048 and effective beam EWHM = 5 
arcmin. 

In particular, each component separation algorithm 
processed differently the multi-contaminated sky regions 
obtaining, as a consequence, different final foreground- 
cleaned regions for each procedure. These regions are 
defined outside the so-called Component Separation Con¬ 
fidence masks (briefly termed CS-masks), which have 
fsky = 0.85,0.96,0.84,0.82, for the SMICA, NILC, SEVEM, 
and Commander maps, respectively. Each CMB map was 
released with its own CS-mask, outside which the corre¬ 
sponding CMB signal is considered statistically robust. 


^ http://planck.mpa-garching.mpg.de/cmb/fnl-simulations/ 

^ Based on observations obtained with Planck 
(http://www.esa.int/Planck), an ESA science mission with 
instruments and contributions directly funded by ESA Member 
States, NASA, and Canada. 


Regarding the masks, it is worth mentioning that the 
Planck team produced the UT78 mask, which is the union 
of the above mentioned masks, a more restrictive one, 
since it has fsky — 0.78, and adopted as the preferred 
mask for analysing the temperature maps [54] . 

All the analyses here are performed upon the products 
of the Planck’s second data release. Moreover, we anal¬ 
yse both the simulated and the four foreground-cleaned 
Planck CMB maps using two masks, namely, UT78 and 
the Commander CS-mask (hereafter. Commander-mask). 
Since the interesting region of the angular power spec¬ 
trum for the current analysis corresponds to ^ = 10 — 50 
(see Eig. [^, we performed some preliminary steps upon 
each Planck CMB map before proceeding. These steps 
are: (1) to degrade the maps to = 512, (2) to cal¬ 

culate the angular power spectrum of each one, and (3) 
to generate the CMB maps from the selected range of 
multipoles. These final maps are the ones we analyse in 
the next section. 



FIG. 3. The multipoles range of the seeded angular power 
spectra for the MC CMB maps, namely, ACDM (solid line). 
Model A (dashed line), and Model B (dot-dashed line). The 
blue, orange, red, and green curves correspond to the spectra 
calculated from the SMICA, SEVEM, NILC and Commander maps, 
respectively. 


IV. RESULTS AND DISCUSSION 

Considering n different thresholds u = ... Un, 

previously defined dividing the range —nrnax to i^max 
in n equal parts, we compute the four MEs {Vk, k = 
0,1,2,3} = (Vb: Fi, V 25 Vs) for the i-th simulated map, 
with i = 1, 2,..., 1000. Then, for the i-th map and for 
k-th ME we have the vector 

^k — (^(^ 1)5 ^^(^ 2)5 •••F^/c(^n))Ifor the z-th map 
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In this work the values chosen for such variables are 
= {3.5,26} (for details, see [50]). 


A. Performance of the MFs as non-Gaussian 
estimators 

Firstly we investigate the performance of the four MFs 
by quantifying how large (in relative difference) is the 
deviation signal in the presence of local-NG in simulated 
maps, as compared with the results obtained in the anal¬ 
ysis of the set of Gaussian ACDM CMB maps. For this 
reason we consider sets of simulated CMB maps with 
the same type of local-NG, but different amounts of such 
contribution, namely with /nl = 38 and /nl = 100, sets 
(iv) and (v), respectively. A similar situation happens for 
the MC sets (ii) and (iii), corresponding to Models A and 
5, respectively. 

We show in Fig.j^our analyses of the four MFs, includ¬ 
ing a comparison using two masks: the UT78 mask (left 
column) and the Commander-mask (right column). From 
top to bottom we analyse the MFs: Area, Perimeter, 
Genus, and ^dusters^ respectively. In all the panels we 
show the mean value of the corresponding MF obtained 
for each of the five sets of 1000 MC maps (we illustrate 
these curves using different colours). We also plot the 
shadow region corresponding to la values of the Gaus¬ 
sian ACDM, i.e., MC set (i). At first sight. Area and 
Genus seems not sensible to the different non-Gaussian 
contents in the sets of MC CMB maps. From other side, 
for the Perimeter and ^dusters MFs, and for the Genus 
case when the Commander-mask is used, some curves ap¬ 
pear close and indistinguishable. For this we find inter¬ 
esting to plot the relative difference between the mean 
curves of the MC sets (ii) to (v) minus the mean curves 
obtained from the Gaussian ACDM maps, as illustrated 
in Fig. 

In Fig. I^the left column refers to analyses with the 
UT78 mask, while the right column refers to analyses with 
the Commander-mask, just to compare the effect left by 
the use of different masks. The four rows of such fig¬ 
ure, from top to bottom, refers to the four MFs: Area, 
Perimeter, Genus, and ^dusters^ respectively. In each 
plot the four coloured curves correspond to the relative 
difference between the mean curves, plotted in Fig. 
of the corresponding MFs from simulated CMB sets (ii), 
(iii), (iv), and (v) minus the Gaussian ACDM MC set (i): 
violet [set (v) - set (i)], green [set (iv) - set (i)], orange [set 
(iii) - set (i)], and red [set (ii) - set (i)], respectively. It is 
worth to notice that the vertical scale is not the same for 
all plots. In the case of the Area-MF, the vertical scale 
is two orders of magnitude less than the other MFs. 

The analyses of these results lead to three important 
conclusions. First, a single statistical estimator is not 
sensible enough to detect any non-Gaussian contribution 
(or a combination of them) present in CMB maps; that 
is to say: various estimators are more efficient to per¬ 
form this task. Second, analyses with different cut-sky 


masks reveal basically the same information, except for 
the Genus-MF as noticed in the plots of the third row of 
Fig. HI Third, and more interestingly, the Gaussian devi¬ 
ations shown in the plots of Fig.j^makes clear that the set 
of MFs are able to discriminate distinct NGs imprinted 
in CMB maps: the non-Gaussian signatures appearing 
in the sets of MC CMB maps emerging from Models A 
and B are not of the type, neither intensity, expected in 
maps contaminated with primordial loeal-NG^ deserving 
a closer analyses (to be done in the next subsection). 


B. Non-Gaussian features in MC CMB maps from 
Models A and B 

Our analyses here are aimed to answer the question: 
are there distinguishable non-Gaussian features left in 
simulated CMB maps generated by Models A and BA 
angular power spectra? For this, we examine the pecu¬ 
liar imprints left in the sets of CMB maps produced from 
the power spectra of Models A and 5, features that can 
be revealed by the MFs in comparison with the Gaussian 
CMB maps generated using the ACDM power spectrum. 
Remember that the power spectra of these Models were 
obtained through the CAMB code assuming a step-like 
potential in the inflation model (these Models differ one 
to the other in the parameter values of the step-like infla¬ 
tionary potential). It is also worth to notice that Model 
A fits the CMB angular power spectrum from Planck 
data better than the ACDM model does (see Fig. 3 in 
Ref. m, and analysis therein); from other side, the an¬ 
gular spectrum from Model B fits the Planck data worse 
than the ACDM model does. 

Our analyses of the four CMB data sets, that is, the 
MC sets (ii) to (v), show consistency at a la level with 
the MFs obtained from the Gaussian ACDM maps, as 
seen in Fig. For this, to go further in the current 
scrutiny we explore the relative difference curves between 
the mean MF from a given MC set, from set (ii) to set 
(v), and the corresponding mean MF obtained for the set 
(i), that is the set of Gaussian ACDM CMB maps. Our 
examinations show, except in the Area-MF case, that 
the CMB maps seeded by the power spectra of Models 
A (set (ii)) and B (set (iii)) evidence a larger Gaussian 
deviation, with fully distinct signature, as compared with 
the results from the MC CMB maps contaminated with 
primordial /oca/-NG, with /nl = 38 (set (iv)) and /nl = 
100 (set (v)). 

In fact, one verifies in the panels of the second, third, 
and fourth rows of Fig. that a distinguished departure 
from Gaussianity comes from the analyses of these CMB 
data produced from Models A and R, with the remarkable 
feature that for the Model B case the Gaussian deviation 
appears with the same signature than for the Model A 
case but it is notoriously larger. In other words, infla¬ 
tionary Models A and B produce a noticeable (in type 
and intensity) Gaussian deviation as measured by the 
MFs, where they are best revealed by the Perimeter and 
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^clusters MFs, and less efficiently by the Genus, beside 
being robust under different masks applied. 

All these results let us to conclude that the non- 
Gaussian signals found in the MG maps seeded by the 
GMB power spectra from inflation Models A and B are 
contributions with a definite signature originated in the 
inflationary phase, suggesting an effective probe to the 
primordial universe. 


C. Mapping non-Gaussian signatures in Planck 
maps 

The third part of our analyses concerns the calcula¬ 
tion of the MFs from the four foreground-cleaned Planck 
GMB maps, considering the above mentioned masks. In 
Figs. 1^ to 1^ we present the relative difference between 
the MFs calculated for each one of these Planck maps 
and the corresponding MF mean obtained for the five 
sets of MG GMB maps. In Tables]^ (for the UT78 mask) 
and [n| (for the Commander-mask) we show the values, 
for 25 degrees of freedom, corresponding to the best-fit 
between the MF from the Planck map and the mean 
MF for each of the five MG sets in scrutiny. First of 
all, our results confirm the point raised above regard¬ 
ing the performance between the four MFs. As a mat¬ 
ter of fact, the Area (Fig. and Genus (Fig. func¬ 
tionals appear to be the less conclusive, each one pre¬ 
senting quite similar differences for all the four Planck 
maps and for both masks. Regarding the Perimeter 
(Fig.[7| and Ndusters (Fig. the panels show a great 
concordance between Gaussian and non-Gaussian (i.e., 
/nl 7 ^ 0) AGDM seeded maps, presenting just tiny rel¬ 
ative differences for the four Planck maps. In these fig¬ 
ures we also observe that comparing the cases regarding 
Models A and 5, they exhibit the same signature but 
very different intensity, being clear that the Perimeter 
and ^clusters functionals from the Model A fits better the 
four Planck maps. Additionally, this last result is corrob¬ 
orated by the analyses, being valid for both masks (see 
Tables [T| and [n|). In other words, our MF’s analyses reveal 
that the non-Gaussian signatures observed in the four 
foreground-cleaned Planck maps are better represented 
by those appearing in the set of MG GMB maps gener¬ 
ated from the Model A angular power spectrum (i.e., the 
MG set (ii)). 


V. CONCLUDING REMARKS 

The main purpose of this work is the analyses of non- 
Gaussian signatures originated in the inflation models 
with a step-like feature in the primordial potential. Such 
models produce oscillatory features in the primordial 
spectrum of scalar perturbations, whose effect on the 
GMB temperature fluctuations can be tested through the 
analyses of simulated GMB maps seeded by the power 
spectrum generated by this model. Here we employed 


TABLE I. values calculated from the MFs (Area (vo), 
Perimeter (vi). Genus (v 2 ), and Ndusters (va)) curves ob¬ 
tained for each Planck map and the average MFs over 1000 
realisations of the sets: (i) to (v). All the MFs considered 
here are calculated using the cut-sky corresponding to the 
UT78 mask. The number of degrees of freedom is 25. 


MF 

Planck 



Data set“ 


map 

(i) 

(ii) 

(iii) 

(iv) 

(v) 


SMICA 

34.0 

35.9 

39.8 

36.7 

42.1 

vo 

NILC 

34.3 

36.3 

40.3 

37.3 

43.0 

SEVEN 

33.8 

35.6 

39.6 

36.0 

41.0 


Commander 

33.0 

34.9 

38.8 

35.5 

40.9 


SMICA 

26.2 

21.0 

78.4 

29.5 

34.1 

Vl 

NILC 

27.1 

22.9 

82.3 

30.5 

35.3 

SEVEN 

24.8 

24.9 

92.3 

27.7 

32.2 


Commander 

27.2 

23.4 

83.5 

30.5 

35.1 


SMICA 

44.6 

43.8 

47.3 

48.4 

51.3 

V2 

NILC 

47.7 

47.8 

50.3 

54.5 

58.2 

SEVEN 

40.4 

39.2 

43.0 

43.5 

45.8 


Commander 

44.7 

43.5 

46.7 

48.0 

50.4 


SMICA 

22.0 

23.5 

57.9 

23.7 

26.3 

V3 

NILC 

27.1 

26.1 

55.3 

29.8 

32.5 

SEVEN 

18.6 

18.4 

50.3 

18.7 

20.1 


Commander 

30.0 

30.9 

63.9 

30.1 

31.7 


“The data sets correspond to: MC CMB maps seeded by (i) 
ACDM, (ii) Model A, and (iii) Model B angular power spectra, 
and two other ACDM seeded sets, but including contributions 
of loeal-NG, namely, (iv) /nl = 38, and (v) /nl = 100. 


four MFs to search for possible Gaussian deviations in 
such synthetic GMB maps. 

Using the five sets of MG GMB maps, and examining 
the effect of using two different cut-sky masks, our anal¬ 
yses initiate by testing the sensitivity of the four MFs to 
loeal-NG, with two different amplitudes, and that NGs 
probably produced by the inflationary Models A and B. 
A first conclusion is that, between the four quantities, the 
Genus-MF is the main affected by the cut-sky. We also 
found that the MFs are able to discriminate the different 
types and amplitudes of NGs. This statement is firstly 
observed from the results presented on Figs. and 
which display the different features of the loeal-NG com¬ 
pared to that ones imprinted on GMB by Models A and 
B. In fact, for the current analyses, we find the following 
order in sensitivity for the MFs: Perimeter > ^clusters 
> Genus Area, from the most to the least sensitive, 
respectively. This result is inferred from the amplitude of 
the corresponding relative differences observed in Fi g. 
as well as from the values shown in Tables and [liP 

The detection of Gaussian deviations, as evidenced by 
departures of the MFs curves compared with the corre¬ 
sponding ones from Gaussian GMB maps, clearly suggest 
the presence of non-Gaussian signals. Exploring this in¬ 
ference more carefully, we analysed the relative differ¬ 
ence between the mean MFs obtained from the four sets 
of maps, sets (ii) to (v), and that one corresponding to 
Gaussian simulations, set (i), seeded by the AGDM power 











TABLE IL values. Same as in Table but for MFs calcu¬ 
lated using the cut-sky corresponding to the Commander-mask. 
The number of degrees of freedom is 25. 


ME 

Planck 



Data set 


map 

(i) 

(ii) 

(iii) 

(iv) 

(v) 


SMICA 

31.9 

33.5 

37.5 

34.0 

38.9 

Vo 

NILC 

32.6 

34.3 

38.4 

35.1 

40.3 

SEVEN 

31.5 

33.1 

37.1 

33.0 

37.3 


Commander 

31.1 

32.6 

36.6 

32.8 

37.3 


SMICA 

28.8 

19.4 

77.3 

30.4 

34.6 

Vl 

NILC 

29.7 

21.5 

82.0 

31.3 

35.8 

SEVEN 

27.7 

22.9 

89.4 

28.8 

32.7 


Commander 

29.9 

22.4 

83.8 

31.4 

35.5 


SMICA 

20.5 

10.1 

18.5 

20.9 

22.5 

V2 

NILC 

35.2 

23.8 

28.1 

37.4 

40.2 

SEVEN 

25.6 

13.6 

19.8 

25.2 

26.3 


Commander 

28.4 

17.3 

23.5 

28.6 

30.2 


SMICA 

13.2 

13.5 

48.2 

13.8 

15.7 

V3 

NILC 

27.3 

23.0 

49.9 

28.6 

31.5 

SEVEN 

13.8 

12.8 

45.6 

14.0 

15.7 


Commander 

18.3 

19.6 

55.2 

18.4 

20.2 


spectrum. From the obtained results, it is inferred that 
MC CMB maps generated from Models A and B power 
spectra present a clear Gaussian deviation, better de¬ 
scribed by the perimeter-MF analyses as sho wn in Fig. 
and quantified by the values in Tables [l| and |liP| whose 

source is not other than the known ‘features’ that char¬ 
acterise these models (see Figs. and [3), therefore orig¬ 
inated in the primordial inflationary phase. Moreover, 
it was possible to reveal the distinctive signature and 
amplitude of the Gaussian deviations present in the MC 
CMB maps seeded by the power spectra from Model A, 
that one that fits Planck data better than ACDM model, 
and Model when compared to the non-Gaussian sim¬ 
ulations (i.e., /nl 7 ^ 0). 

Finally, in the last part of this work we analysed the 
four foreground-cleaned Planck CMB maps in the current 
context. The analyses of the relative differences between 
the MFs of each Planck map and the average MFs from 
the five sets of MC CMB maps confirm a good agree¬ 
ment between the Planck data and data coming from the 
MC seeded by Model A power spectrum (i.e., the MC 
set (ii)). In addition, we use the goodness-of-fit test to 
get an overall measure of this concordance obtaining the 
values presented in Tables and which confirm 
this result. This lead us to our main conclusion con¬ 
cerning Planck data: non-Gaussian signatures present 
in the four foreground-cleaned Planck maps are better 
represented by the MCs CMB maps seeded by Model A 
power spectrum than is done by those MCs seeded by the 


^ Notice that the MFs are not equally sensible to all types of NG 
(as illustrated, for instance, in the panels of Fig. [^. Moreover, 
the MFs look for all types of NGs, primordial and secondary 
ones, present in the CMB maps in analyses, therefore one does 


ACDM (Gaussian and loeal non-Gaussian simulations) 
and Model B power spectra. 

Considering that any convincing deviation from Gaus- 
sianity results from a robust detection by the MFs, which 
include analyses with different masks. In addition, such 
departures are larger when the power spectrum features 
appears exaggerated (this refers to Model B ‘features’ 
with respect to those of Model A, see Figs. and [^, 
then the main conclusions of the current analyses regard¬ 
ing the MC CMB maps derived from inflationary models 
with a step-like potential (i.e., the MC sets (ii) and (hi)) 
are: 

• The non-Gaussian contributions found in these sets 
of MC CMB maps are of primordial origin. 

• These net Gaussian deviations, as robustly de¬ 
tected by the MFs, are not accounted by primordial 
/(9ca/-NG, neither in intensity nor in signature. 

• The non-Gaussian signatures in Planck CMB data, 
as revealed by MFs, are better described by MC 
maps seeded by ModeVs A CMB spectrum. 

Additionally, it is opportune to inquire about the type 
and level of the NG detected. According to the results 
of Chen et al. [23], inflationary models with a step-like 
potential would produce detectable levels of NG of equi¬ 
lateral type, with expected amplitude — 10, for 

the step parameters considered in Model A. Moreover, 
as showed in Table [l| the best-fit satisfied by the 
MCs from Model A indicates that the NG observed in 
the foreground-cleaned Planck maps is fully compatible 
with the NG (of equilateral type as predicted by [23] 1 
found by our perimeter-MF analyses in these MC CMB 
maps (see Fig.[^. These facts, altogether, show the good 
agreement between the equilateral NG measured by the 
Planck collaboration, namely, = —16 ±70 [21], and 

the NG predicted by [23] for Model A and detected by 
the perimeter-MF in the four Planck CMB maps. 
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FIG. 4. Comparison between the average Area (vo), Perimeter (vi), Genus (V 2 ) and, Ndusters (vs) MF vectors over 1000 
realisations of each dataset, using the two masks, UT78 (left) and Commander-mask (right). The blue, red, and orange curves 
corresponds to the average MF obtained from dataset seeded by the ACDM, Model A and Model B power spectra, respectively. 
The green and violet curves represent the average MFs calculated from MG CMB maps seeded by ACDM power spectra and 
contaminated by loeal-NG with /nl = 38 and 100, respectively. The gray region represents the la level calculated from the 
MFs curves obtained from the ACDM seeded dataset. 
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FIG. 5. Relative difference between the mean curves plotted in Fig. of the corresponding MFs for the sets of CMB maps 
(ii) to (v) minus the mean curves obtained for the Gaussian AGDM set (i) (see the text for details about these hve sets of 
MG GMB maps). The left (right) column refers to analyses with the UT78 mask (Commander-mask), to compare the effect left 
by the use of different masks. From top to bottom the four rows of panels refers to the four MFs: Area, Perimeter, Genus, 
and ^clusters, respectively. In each plot the four coloured curves correspond to the difference between the MFs from simulated 
GMB sets (ii), (hi), (iv), and (v) minus the AGDM set, i.e., the Gaussian MG set (i): violet [set (v) - set (i)], green [set (iv) - 
set (i)], orange [set (hi) - set (i)], and red [set (ii) -set (i)], respectively. It is worth to notice that the vertical scale is not the 
same for all plots. In the case of the Area-MF, the vertical scale is two orders of magnitude less than the other MFs. 
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FIG. 6. Relative difference between the Area-MF vectors (vo) calculated from each Planck map (from top to bottom: SMICA, 
NILC, SEVEN and Commander maps, respectively) and the average Area-MF ((vo)) obtained from the five sets of MC CMB maps, 
namely, ACDM seeded Gaussian (blue) and non-Gaussian (/nl = 38 and 100; green and violet, respectively) maps, and Model 
A (red) and Model B (orange) seeded maps. The left and right panels correspond, respectively, to the analyses performed using 
the masks UT78 and Commander-mask. 
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FIG. 7. Same as in Fig. but for the Perimeter-MF vectors. 
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FIG. 8. Same as in Fig. but for the Genus-MF vectors, 
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FIG. 9. Same as in Fig. but for the Ndusters-MF vectors, 


















